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Abstract 

In this paper, we study the Ginzburg-Landau equations for a two 
dimensional domain which has small size. We prove that if the domain 
is small, then the solution has no zero, that is no vortex. More precisely, 
we show that the order parameter 4" is almost constant. Additionnally, 
we obtain that if the domain is a disc of small radius, then any non 
normal solution is symmetric and unique. Then, in the case of a slab, 
that is a one dimensional domain, we use the same method to derive 
that solutions are symmetric. The proofs use a priori estimates and 
the Poincare inequality. 

1 Introduction and main results 

In this paper, we study the properties of a superconducting cylinder sub- 
mitted to an exterior magnetic field Ho parallel to the axis of the cylinder. 
According to the Ginzburg-Landau theory of superconductivity, the sample 
is in a state that minimizes the following energy: 

4(«,A)= / |(iv-iA)*| 2 + i(|*| 2 -l) 2 + (curl A -H ) 2 dfi. (1.1) 

Here f2 is a simply connected domain in IR 2 with characteristic length d. 
As usual, k is the Ginzburg-Landau parameter that characterizes the type 
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of superconductor, A(x, y) is the vector potential, so that curl A is the 
magnetic field and ^f(x,y) is the order parameter. Because Ho = Hqg z 
is along the z axis, we can assume without loss of generality that A = 
(A\ , Ai , 0) . For a detailed description of the model, one may refer to Q or 

0- 

The minimization process yields the classical Ginzburg-Landau equa- 
tions (see || for more details): 

(iV-fA) 2 ^ = ^(|^| 2 - 1) in O, 



-curl curl A = ^(**V$ - Wf*) + A|f| 2 in fi 
|| - MA • n = on 90, 
curl A x n = Hq x n on <9$7. 



It is common to use the gauge where div A = and A-n = on <9f2. Notice 
that this system has a special set of solutions with ^ = and curl A = Ho, 
called normal solutions. They correspond to a case where superconductivity 
is destroyed. In the following, we will be concerned with solutions which are 
not in this set. 

In this problem, there are 3 variable parameters: k, Hq and d the charac- 
teristic size of Vl. According to the values of these parameters, the properties 
of solutions (existence of non normal solutions, number, symmetry) change. 
It has been the purpose of various authors to study the properties of solu- 



tions in some asymptotic limits: for instance, Giorgi and Phillips [10] have 
proved that for £1 and k fixed and Hq large enough, then the only solution is 
the normal solution. In the case where d is fixed and k is large, Sandier and 
Serfaty [|D|] have proved that the solutions can have a lot of vortices for a 
certain range of magnetic field. There have been many other contributions 
in other asymptotic regimes. 

In this paper, we are interested in the case where d, the characteristic 
length of Q is small, in particular compared to 1/k. In this setting, it is 
expected that \& is almost constant and has no zero, that is no vortices, since 
a vortex is of size 1/k. This result is often used in the physics literature, 
but we did not know any rigorous proof. This is what we show: 

Theorem 1 Assume that D is a fixed simply connected bounded domain 
and let SI = dD. For any do > 0, there exists d% > 0, such that if d < 
min(do, di/n), for any Hq = HQ(d, k), then any solution of(GL) which is not 
a normal solution is such that has no zero. Moreover, |\E , /| , I / | 00 — 1| < Cud. 

We let Hq be a function of d and k because it is expected that a supercon- 
ducting solution exists up to fields Hq of the order of C/d as we will explain 
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in more details later. Additionally, we prove a symmetry and uniqueness 
result when the domain is a small disc: 

Theorem 2 Let Q be a disc of radius d. There exist constants do and d\, 
such that if d < min(do, g?i/k), for any Hq = Ho{d,n), then any solution of 
{GL) is radially symmetric, that is ^{x,y) = \I/(r) and A{x,y) = A(r)eg, 
where r = \J x 2 + y 2 . 

Theorem 3 Under the hypotheses of Theorem |], if there exists a non nor- 
mal solution (vP, A) of {GL), then this solution is unique up to multiplication 
of ^ by a constant of modulus 1. 



Corollary 4 Under the hypotheses of Theorem^, there exists H* = H*{d, k), 
with lim^_»o dH*{d, k) = 2\/2 and for Hq < H*, there exists a unique non 
normal solution of {GL), while for Hq > H*, the only solutions are the 
normal solutions. 

We make a change of variables x' = x/d, y' = y/d so that the new 
variable lies in a domain of unit size D. We also define B = KdA and 
ho = nd 2 Ho- In the following, we will assume additionnally that 

div B = in ft and B n = on d£l. (1.2) 

The equations then become 

(V -iB) 2 f = K 2 d 2 ^j{\^\ 2 -1) in D, 
AB = d 2 (|(V*VV- W")+B|^| 2 ) in D, 
f = on dD, {GLd) 
curl B x n = ho x n on dD. 

Note that another way of writing the equation for B is 

AB = —d 2 {iip,Vip — iBip), (1.3) 

where ( . , . ) is the real part of the scalar product in C. 

We allow ho to vary with d and k but we will prove that if there exists 
a solution, then in fact ho is bounded, that is Ho is bounded by C / ad 2 . 

The proof consists in obtaining a priori estimates for the solutions (ip, B). 
This is done in Section 2. In Section 3, we use these a priori estimates and 
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the Poincare inequality to derive that ip is nearly constant, hence has no zero. 
This will prove Theorem |]. Then, in Section 4, we define the functions 

$(x, y) = ip(-x, y), B(x, y) = 

w(x, y) = — — (i>{x, y) — Tp(x,y)), z(x) = B(x, y) - B(x, y), (1.4) 

\V\oo 

which satisfy elliptic PDE's with small right hand side terms. Then we use 
that ip is nearly constant and hence we get that w and z are identically 
zero. This will prove Theorem ^. In Section 5, we obtain the uniqueness 
result proving that any solution is necessary a local minimizer of the energy. 
Finally in Section 6, we show how our proof can be adapted to the one 
dimensional case to provide a symmetry result in this setting and in Section 7 
we summarize our results in terms of bifurcation curves and prove Corollary 
I 

Remark: If 0, is simply connected and symmetric in the y direction, then 
our proof also gives that ip and B are symmetric, in the sense that ip = ip 
and B = B. 



-B x (-x,y) 
B2{~x,y) 



We note that in the case with no magnetic field, Bethuel, Brezis and 
Helein Q have proved the symmetry of solutions when the domain is small 
using the Poincare inequality. The radial symmetry of local minimizers has 
also been studied in [^] for a different system with no magnetic field. Here 
the proof is more involved since we have to deal with the system with applied 
magnetic field. 

Let us point out that the uniqueness proof of Section 5 is inspired by || 
where it is obtained that any radial solution with a giant vortex of degree 
iV at the origin is a local minimizer of the energy. 



2 A priori estimates 

Proposition 5 Fix p > 1. Assume that (psi,B) is a non normal solution 



of (GLd) and (lji). For all constants do and d\, if d < min(cio, di/n), 
then ip and B are bounded in W 2,P (D) by constants independent of d and k. 
Moreover, for fixed k, 

lim ho = and lim ||B||i / i/2,pi'n , i = 0. (2.1) 
The goal of this section consists in proving Proposition [f| 
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2.1 First estimates for ip 

First recall from Du-Gunzburger-Peterson @ that 

< 1 a.e. (2.2) 

Next we define 

n(x,y) = \il>{x,y)\ 2 and II = V-iB. (2.3) 
One can compute easily that 

^Au + K 2 d 2 u(l -u) = |IIV>| 2 . (2.4) 

Integrate fl2.4j ) in D, use the boundary condition dip/dn = and get the 
key estimate for ip: 

\\UiP\\ LHD) < KdM\^ {D) . (2.5) 

2.2 Estimates for B 

We decompose B = Ziobo + b r where bo is chosen such that 

curl curl bo = and div bo = in D, curl bo = e z and bo • n = on dD, 

and b r = B — /lobo- 

A particular choice for bo when $7 is a disc is bo = |(— y,x,0). Note 
that in any case bo and b r are vectors in the x — y plane and their curl are 
in the z direction. Making the difference between the equations for B and 
bo, we obtain the equation for b r 

curl curl b r = d 2 {iip, II^) and div b r = in D, (2.6) 
curl b r = and b r • n = on dD. (2.7) 

Since div b r = and b r ■ n = on dD, then ||curl b r ||£ P (£)) is a norm in 
W 1 ' P (D), that is 

l|br||^i,p(-D) < C||curl b r || LP(D ). (2.8) 
Now Cauchy-Schwarz inequality and (|2.6[) imply that 



|curl curl b r | 2 < d 2 \Uip\ 2 . (2.9) 



Note that | . | 2 is the 2- norm for vectors in K 2 . From fl2.9| ), we get the W 2 ' 2 
bound for b r . Thus, because of (^), ( |2.9| ) and the boundary condition 
( |2.7D , it follows from the Poincare inequality that 

1 1 curl h r \\ w i,2( D ) < Cud 3 . (2.10) 
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We gather (|2.8[ ) and ( 2.1 0| ) to obtain the key estimate for b r : 



||b r || W 2,2 (D ) < Cnd 6 , (2.11) 

where C is independent of d and k. In order to get the bounds for B, we 
only need to prove that lim^o ho = 0. Let us assume that ho is bounded 
below by some constant m. Then we use an estimate by Giorgi and Phillips 
JiCj| (Lemma 2.8 p. 349): there exists C2 (which depends on m) such that 

C 2 ho [ |Vf</ |VV-i/i boVf • (2.12) 
Jd Jd 

In order to bound the right hand side of ( [2.12D , we write hoho = B — b r and 
use Q and flUl]) to get 

(C 2 h ) 1/2 ([ |V| 2 ) 1/2 <^(/ |Vf) 1/2 (l + dd 2 ). (2.13) 
Jd Jd 

So we obtain a contradiction with the lower bound for ho when d is small 
enough. Hence B is bounded in W 2 ' 2 , independently of d and k. 

Once we have bounded B, the equation for if; can be thought of as a 
linear elliptic equation with bounded coefficients, so that the classical elliptic 
estimates imply bounds for ip. 

Finally, in order to get the W 2,p bound for B, that is for b r , we bootstrap 
equation (^) (see fH). Thus we obtain that 

1 1 curl b r I |w-2.*>(£>) < Cd 2 . (2.14) 



and (2.1) is true for B. 



3 i/j is nearly constant 

Proposition 6 For all do > 0, there exists d\ > 0, such that if d < 
mm(do,di/K), then ||V'0/|'i/ ; |oo||oo and IV'/IV'loo — 1| are small, thus i\) has 
no zero. 

We define <p = t/VIV'Iocm then the equation for <j) can be written as 
Acj) - 2iB ■ = « 2 dV(MLM 2 - !) + l B | 2( ^ 

and for p < 2 

|B • V4>\ p < |B • n^| p + |B 2 ^»| p < I—,— n^| 2 |B|^p_ + |B 2 | p . (3.1) 

V OO 2-p 
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Since \(j)\ < 1, and because of (2.5) and the previous estimates for B, we 



see that the right hand side of (3T) is small when d is small. Hence we 



go back to the equation for (j) and the Agmon Doughs Nirenberg estimates 
imply a bound for <fi in W 2,p for any p < 2. Hence by Sobolev embedding, 

is bounded in L p for any p and going back to the equation for <p, we 
derive a bound for <f> in W 2,p for any p, hence in C . In particular, V(p 
is equicontinuous. Now, when d tends to 0, we know that HB^ tends to 
0, so if we multiply the equation for <f> by <p* and integrate, we find that 
|| ^7((>\\l 2 (D) is small. Thus, since V</> is equicontinuous, we see that ||V</>||oo 
tends to with d, which proves the first part of Proposition ||. 

In fact we can get a precise estimate for the smallness of (V'/lV'loo — 1|- 
Take the equation (|2.4[), multiply by u and integrate to get 

Vu| 2 < K 2 d 2 j u 2 {l - u). (3.2) 

D JD 



Next we define v = u/lip] 2 ^. Then (|3.2|) implies that 

\Vv\ 2 < K 2 d 2 \D\. 



D 

We call 

v mean the mean value of v in D. The previous estimate and the 
Poincare inequality yield 

\v ~ V m ean\\L 2 (D) < C Kd. (3.3) 

Let us call r an upper bound for Vi>, independent of d and k. We have seen 
in the first part of the proof that it exists. 
Let us now prove that 

\v(x) — v mean \ < a \fx G D, where a 2 = lQCornd/^fn, (3-4) 

where Cq comes from fl3.3|) . Assume that ( |3.4 ) is not true. Then there exists 
a point xq where this does not hold. Set r/ = a/Ar. Note that when nd is 
small, then r] is small. If dist(xo, dD) < 2r/, then for nd small enough, there 
always exists a point x\ in D with dist{x\,dD) > 2rj and dist(xi,a;o) < 2rj. 
If dist(xo, dD) > 2rj then we set x\ = xq. Now we have \v(x) — v mean \ > a/4 
for \x — x\\ < T). Define D% = D(xi,rj), then 

\\v - v mean \\ L 2( Dl ) > rfy/na/A. 

This provides a contradiction with the definition of a and ( |3.3|) . 

Now recall that v is necessarily 1 for at least one point in D. So we use 
Q to get 

|^0r)/|V>|oo - 1| < Cnd Vx £ D. 
This implies in particular that for small nd, ip is never equal to zero. 
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4 Radial symmetry of the solutions in the case of 
a disc 



In this Section, we assume that f2 is a disc of radius d. We are going to use 
the previous a priori estimates to obtain radial symmetry. We define the 
functions w and z as in (|1.4| ). They satisfy the following equations 

2i 

Aw - 2iB ■ Vw - T-^z • Vip = K 2 d 2 w(\il)\ 2 + M + - 1) 

moo 

+w\'B\ 2 + 7-7— z • (B + B) in D, (4.1) 



^ = on 3D, (4.2) 
on 

curl curl z = d 2 [(i|V>|ooW) n^) + (iip, \ip\ooVw — iTS\ip\ooW — itpz)] 

and div z = in D, (4.3) 
z • n = and curl z x n = on 3D. (4-4) 



Proposition 7 There are constants do and d\ such that, if d < min(do, 
di/n), then w = z = 0. 

Let us multiply (4.1) by w and integrate. The first term only gives an 
integral of the gradient of w because of the boundary condition (4.2). We 
recall that < 1 and \w\ < 2 to obtain 

/ \Vw\ 2 < 4||B||oo / \Vw\ + {4 K 2 d 2 + ||B||^) / w 2 

JD JD JD 

+4(1^-^1100 + HBIU) f 



|z|. 
id 



Then we use our previous estimates for ip and B in Proposition |5| and |(| In 
particular we recall that HBHoo and || p^j— VV>||oo are small when d is small 
to get: for all e > 0, if d is small enough, then 

[|Vm>[|2 < £ IM|2 -I— e" 1 1 z 1 1 2 - 

Because of the definition of w, one can see easily that w has mean value 
zero. We can use the Poincare inequality to get 

AilMb < ^IMk + e||z||2, (4.5) 



S 



where \i > is independent of d and k. Thus 

e 

i» 2 < z 2 ^ z 2, (4.6) 

/i — £ 

provided s/(fi — e) is less then 1. 

Similarly, we multiply equation (4.3) by z, integrate and use that the 
boundary condition implies curl z = on dD to get 

/ | curl z\l < d 2 / |z| 2 |iT/> + Vw + B + z\ 2 }. 
Jd Jd 

We can apply the previous estimates which give bounds on ip and B to 
derive 

||Vz|| 2 < e||z|| 2 . (4.7) 

Then we write z = z m + zt where z m is the mean value of z (so z m is a 
vector) and zj has mean value zero. Again we use the Poincare inequality 
for zt and equation (|4.7|) to get 



£ £U 

% 2 < z m 2 and \\w\\ 2 < -, z m 2 - (4.8) 

H - £ {n - e) 2 



Now we integrate ( |4.3| ) over D, use the boundary conditions which imply 
that the curl curl term contribution vanishes and we get 



f -i^Lz= / (iw,nip) + (iij),Vw-iBw). (4.9) 
Jd \tp\io Jd 



This yields 



\ z m\ [ S- < Hl2||n^|| 2 + W2||V«;-iBi(7|| 2 + / \z t \<Ce\z m \. 
Jd mio Jd 

The last inequality comes from ((T^) and the bounds for ip and B in Section 
3. If z m is not zero, the boundedness of Hip, B and (|4.8[) imply that 



D 



is small. This provides a contradiction with Proposition ||, where we have 
proved that IV'l/IV'loo is nearly 1 for small d. 

Let us complete the proof of Theorem ||. Indeed, once we know that w 
and z are identically zero, then it implies as in M that if) and |B| are radial, 
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where | . | is the modulus. Thus, we see that either B vanishes everywhere 
on the circle |(x,y)| = r and there is nothing to prove, or nowhere on 
this circle. In this case, we can argue that B(r, 9) is in the direction of 
eg = (— s'mO, cos 9). Indeed, since z = 0, it implies in particular that 
Bi(r, 0) = 0, hence B(r, 0) is in the direction of (0,1). Then we use a 
reflection in the axis T in the direction of (cos 9, sin 9) and derive that the 
corresponding w and z are zero. This implies that B(r, 9) ■ e r = 0. Since 
B\(r, 0) = 0, a continuity argument implies that B(r, 9) = B(r)eg, as we 
had claimed. 

Let us show that one can choose ip to be real. Since ip is radial and 
has no zero, we can write tj) = f{r)^^ r \ where / is real. The estimate of 
Proposition ^ yields that 4>'{r) is small for small d. Thus, if we change gauge 
from (i/>,B) to (/, Q) where Q = B — (p'(r)e r , then (/, Q) is a solution of 
(GL d ), Q • n = on dVL and |Q| oo is small. Notice that in this section, we 
have not explicitely used that div B = 0. We have only used that |B|oo is 
small, which was proved in Proposition || under the hypothesis div B = 0. 
Now, since |Q|oo is small, this is enough to do the proof of this section with 
(/, Q) instead of (?/>, B). It yields that Q is along eg and that eft' = 0. Thus 
one can choose tp up to multiplication by a constant of modulus 1. This 
finishes the proof of Theorem ||. 

5 Uniqueness for the ball 

In this section, we assume again that D is a ball. We have proved in the 
previous section that we can assume that ip = f(r) is real and B = Q = 
Q(r)eg. Then the Ginzburg-Landau energy in this case is 

E(f, Q)= I f + fQ 2 + ^f(f- l) 2 + -kW - h f. (5.10) 
Jd 2 d z r 

We also have div Q = then ||^(^(5) / ||l2 is an H 1 norm. We are going to 
prove that any radial solution of (GLd) is a local minimizer of the energy 
( |5.10| ). If (/, Q) is such a solution, then for any g and P = P(r)eg in H 1 , 

= E'(f,Q).(g,P)= [ f'g' + QPf + Qlfg 

Jd 

+ h\rP)i){\rQ)> -ho) + K 2 d 2 (f -l)fg. (5.11) 

d A r r 

Moreover 

E"(f,Q)-(g,P) = [ g> 2 +g 2Q2 + 4fgQP + pif 
Jd 
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+ ^(-(rP)') 2 + ^W(3/ 2 -l)/5 (5.12) 



Proposition 8 Let (/, Q) be a radial solution of (GL^) with d and dn small. 
Then E"(f,Q) ■ (g,P) > C(\\g/ f\\ H i + ||P||#i), where C depends on d, nd 
and ||/]|oo- 



The proof follows ideas from |J. First compute E'(f,Q) ■ (g 2 /f,0)- This 
yields 



o = / o /'(^)VQ 2 + A 2 (/ 2 -i)s 2 

Hence 

£"(/,Q>( 5 ,P) = / D ((|) , ) 2 / 2 +4W++ J P 2 / 2 +^(^(rP) / ) 2 +2^ ( iV/ 2 . 
We write 

4/ 5 QP + K 2 d 2 g 2 f 2 = (Kdfg + -^ p Q) 2 ~ -^2 p2 Q 2 - 

Hence 

p"(/, o) • G* P) > / D ((f )') 2 / 2 + / 2 + ^(>^)') 2 - ^ 2 Q 2 - 

Now we use that for d small enough, ||Q||oo is small hence 8 J D P 2 Q 2 < 
k 2 J D (^(rP)') 2 , thus the estimate of the Proposition holds. 

This Proposition yields that any solution of {GL^) is non degenerate. We 
are going to prove uniqueness using the bifurcation curve. More precisely, 
any radial solution of (GL^) solves 

with the boundary conditions 

/'(0) = f'(l) = and q(0) = 0, (<?' + ±q) (1) = h . (5.14) 



We use shooting techniques, that is solve ( |5.13 ) with the initial conditions 



/(0) = 0, /'(0) = 0, g(0) = 0, q'(0) = a. (5.15) 
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Then the same kind of proof as in Kwong's paper []IT| which applies the 
Sturm comparison principle to equation ( |5,13|) shows that f(r; (3, a), q(r; (3, a) 
and q'(r;/3,a) are increasing functions of (3 and a. Moreover, for each (3 in 
(0, 1), there exists a unique a such that /'(l; /?, a) = 0, that is the boundary 
condition is satisfied by /. Then 

h((3)=(q' + ±q)(l;P,a((3)) 

is a continuous function of (3 with limit as (3 tends to 1 and a finite limit 
as [3 tends to 0. The nondegeneracy result of Proposition || implies that this 
curve is decreasing and it provides uniqueness. 

Note that uniqueness could also be proved by combining the compact- 
ness of solutions and Proposition || with Crandall Rabinowitz bifurcation 
Theorem. The Crandall Rabinowitz bifurcation Theorem ensures that only 
one branch comes off the normal solutions. Then the non degeneracy result 
of Proposition || and the implicit function Theorem guarantee that the value 
P continues locally uniquely as ho is varied. This ensures uniqueness. 



6 The 1-dimensional case 

When the superconducting material is an infinite slab of thickness 2d be- 
tween the planes x = —d and x = d, it is usual to assume that both \E r and 
A are uniform in the y and z direction, and that the exterior magnetic field 
is tangential to the slab, that is Ho =(0,0, Hq). A suitable gauge can then 
be chosen so that ^ = f(x) is a real function, and A = q(x)e y , where e y is 
the unit vector along the y direction (see |7| for more details). The model 
can then be simplified to a system of 2 coupled ODE's (/, q) satisfy 

^f" = f(f + q 2 -l) m (~d,d), 
f'(±d) = 0, 

q" = qf 2 in (-d,d), 
q'{±d) = H . 



In this setting, there are two types of solutions: symmetric solutions where 
/ is even and q is odd, and asymmetric solutions. A complete numerical 
study of the number and symmetry of these solutions has been done in 
Using the same techniques as in the previous sections, we can prove 

Theorem 9 There exist constants do and d±, such that if d < min(cio, 
di/n), then any solution of (gld) is symmetric, that is f is even and q is 
odd. 
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This theorem together with the result of uniqueness for symmetric solu- 
tions proved in [|j give a global uniqueness result for the solutions of (gld) 
with small d. 

The proof, as in the previous sections consists first in deriving a priori 
estimates for the functions / and q. We recall from Q that a solution is such 
that / has a unique maximum, which we call xo, with = /(xo) = ll/lloo, 
and q is increasing with a unique zero x\. A similar proof to what we did 
in Section 2 and 3 yields: 

Proposition 10 There exist constants do and d\, such that if d < min 
(do, d\j ' k), then 



/3 



1| < C«d 5 and |g(x) - ff ^| < Cd 2 . (6.1) 



Using similar techniques as in one can get more precise estimates, 
which in particular give a relation between (5 and /do. 

Proposition 11 Let < d < min(l, 1/2k), then 

/3(l-y(x-x ) 2 )</(x)</3 VxG(-d,d), (6.2) 

and i/iere exists a such that 

a(x — xi) < q(x) < a(x — xi)(l + 2j3d) Vx G (xi, d), (6-3) 
a(x — xi)(l — 2j3d) < q(x) < a(x — xi) Vx G (— d, xi), (6.4) 

|xi| < 2K 2 d 3 + 2d 2 , (6.5) 

These estimates mean that / is nearly constant and that q is nearly equal 
to ax, because its zero xi is very small and recall from Proposition fill that 
a and Hq are very close. We will not give the proof of Proposition |ll] since 
Proposition is enough for our purposes. 

The second step of the proof consists in defining the functions w and z 
as before (note that the distances have been non dimensionnalized by d as 
in the previous sections: 

w (x) = _L_(/(a.) _ f(-x)) and z{x) = + , ( 6 . 7 ) 
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We derive the equations satisfied by w and z, multiply them respectively by 
w and z and use the estimates of Proposition |l0| to obtain 



\w"\ < e\w\ + e\z\ and \z"\ < e\w\ + s\z\ 



Then we use the Poincare inequality to get the equivalent of (4.6) and (4,7). 
Finally, we see that since /_., f 2 q = 0, then 



1 f 



1 (/ + /) 

qW R~ ■ 

-1 p 



(6.8) 



As previously, we derive that if z m is different from zero, ( |6.§| ) provides a 
contradiction with (6.1), which means that f(x)/f3 is nearly constant. 



7 Remarks on the bifurcation curve 

In order to describe the solutions of the Ginzburg-Landau equations, it is 
common to draw the bifurcation diagram. The natural bifurcation diagram, 
which appears in the physics literature, is to plot the infinity norm of W 
(that we call (3) against Hq. In our setting, recall that is nearly constant. 

For the 1-dimensionnal model, Kwong jll| has proved that for each f3, 
there exists a unique Hq((3) such that (gld) has a symmetric solution. Recall 
that in our setting, we have proved that solutions are indeed symmetric. 
Aftalion and Troy [|J have established a complete picture of the bifurcation 
curves in the different regimes of k and d. In the regime where d is small 
compared to 1/k, then the curve (3 against Hq is decreasing from 1 to 
as Hq increases from to C/d. In Q, it is proved that lim^o dHo((3) = 
\/3(l — (3 2 ) uniformly with respect to (3 in (0, 1). Thus, for Hq larger than 
v3/cZ, the only solution is the normal solution. For Hq smaller, there exists a 
unique solution which is not normal; it is symmetric and has \& or / almost 
constant to [3 between and 1. The value of Hq for which the normal 
solution loses its stability has been studied in detail by Bolley-Helffer || . 

In the 2 dimensionnal case, the same kind of behaviour holds. A natural 
question is to wonder when the superconducting solution exists. Using the 
result of Giorgi-Phillips, we have proved that if there exists a non normal 
solution then d 2 HQ tends to with d. Hence, if Hq is not small enough, the 
only solution is the normal solution. On the other hand, it is known that 
for Hq small, there exists a superconducting solution close to 1. 

When there is a superconducting solution, we have also proved that these 
solutions are such that ^ is nearly constant. Let us call this constant (3. 
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Then there is a relation between the value of j3 and Hq. This relation is 
obtained by integrating the equation for ^ using the boundary condition. 
One gets 

J ^O 2 + \A\ 2 - 1) = 0. 
In the 1-d case, since q(x) is equivalent to Hqx, we get 

lim d 2 Hl = 3(1 -P 2 ), 

d^O 

as proved in 0|. In the radial case, one can prove with a similar estimate to 
Section 2 that A(r) is asymptotically Hqt/2, hence 

limd 2 ^ = 8(l-/3 2 ). 
d^o u 

This gives an asymptotic limit of the bifurcation curve for small d and 
finishes the proof of Corollary ||. 
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